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1  Introduction 

Thermal  neutrons  with  energies  on  the  order  of  0.025eV  are  of  considerable  interest  to 
the  Department  of  Defense  and  for  commercial  applications.  Unlike  high-energy  photons, 
neutrons  easily  penetrate  high  density  targets,  but  get  effectively  absorbed  by  low  density 
materials  like  paraffin,  nylon  or  explosives.  This  makes  them  attractive  complements  to 
X-rays  for  radiographic  applications,  e.g.  for  the  detection  or  inspection  of  explosives 
inside  steel  casings.  The  key  challenge  is  to  develop  a  compact  generator  for  thermal 
neutrons  with  a  flux  of  at  least  10  neutrons/cm  s  which  is  sufficient  for  radiographic 
applications. 

The  limited  available  amount  of  radio-isotopes  like  Californium  252,  combined  with  the 
relatively  short  half-life,  safety  constraints  and  regulatory  requirements  make  them 
unattractive  for  wide-spread  use. 

An  alternative  design  exploits  the  Deuterium-Tritium  (D-T)  fusion,  which  generates 
Alpha  particles  and  fast  neutrons.  In  these  sources,  Deuterium  ions  are  accelerated  to 
about  130  keV  and  hit  a  Tritium  target.  In  order  to  be  attractive  for  radiographic 
applications,  the  14.1  MeV  neutrons  have  to  be  thermalized  in  an  external  moderator,  for 
example  a  paraffin  layer.  This  reduces  the  useful  neutron  flux  by  about  two  orders  of 
magnitude.  D-D  reactions  produce  neutrons  with  lower  energy  (2.4  MeV),  but  one  must 
also  account  for  the  branching  ratio  into  charged  products;  i.e.,  d(d,p)t,  and  the  reduced 
cross-section  compared  to  d(d,n)a. 

The  acceleration  of  Deuterium  ions  is  usually  accomplished  in  a  diode  configuration. 
Recently,  the  interaction  of  high-power  laser  pulses  with  plasmas  or  solids  has  attracted 
considerable  interest  for  the  acceleration  of  electrons  or  ions  [1,  2],  Especially  the 
acceleration  of  ions  to  MeV  energies  opens  entirely  new  possibilities  for  small  scale  ion 
accelerators  and  possibly  for  neutron  sources. 

In  this  project,  we  therefore  investigated  different  target  configurations  for  the  potential 
use  as  neutron  generators  via  numerical  simulations. 

2  The  Plasma  Simulation  Framework  VORPAL 

The  versatile  plasma  simulation  code  VORPAL  [3]  was  originally  designed  for  the 
investigation  of  laser  wakefield  acceleration  of  electrons.  However,  the  code  design  is 
sufficiently  general  that  it  is  capable  to  address  entirely  new  application  areas,  like 
modeling  of  superconducting  RF  cavities,  breakdown  phenomena  in  microwaves  guides 
or  investigations  of  astrophysical  plasmas.  Figure  1  shows  examples  of  such  simulations. 
The  code  has  recently  been  enhanced  by  a  variety  of  models,  including  ionization 
processes. 

The  model  for  ionization  processes  is  a  hybrid  between  classical  particle-in-cell  [4]  and 
the  Direct- Simulation  Monte-Carlo  model  [5],  Particles  are  pushed  in  their  self-consistent 
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electromagnetic  field.  At  every  time-step,  particles  within  a  cell  are  considered  for 
possible  collision  (and  therefore  reaction)  by  Monte-Carlo  methods.  The  main  parameters 
required  in  this  model  are  energy  dependent  cross  sections  for  ionization.  While  the 
model  was  mainly  designed  for  ionization  processes,  it  can  relatively  easily  be  adapted  to 
other  reaction  types,  including  fusion. 


Figure  1:  Left  panel:  Example  of  a  VORPAL  simulation,  (superconducting  RF  cavity).  Right  Panel:  Mono- 
energetic  electron  beam  generated  by  laser-wakefield  acceleration  (Simulation  by  John  R.  Cary, 
Visualization:  Peter  Messmer). 


2.1  Ionization  Model 

The  ionization  model  in  VORPAL  was  successfully  used  to  reproduce  ion  spectra  in 
Electron-Cyclotron  Resonance  Sources.  Figure  2  shows  a  comparison  of  the  kinetic 
impact  ionization  model  in  VORPAL  with  a  spectrum  for  a  hot  population  of  oxygen 
measured  at  VENUS  [7].  For  experimental  reasons,  the  spectrum  was  not  measured  for 
pure  oxygen,  but  rather  a  mixture  of  Oxygen  and  Helium.  The  resulting  spectrum  at 
mass/charge  ratios  of  4  was  therefore  a  blend  between  the  oxygen  and  the  helium 
contribution.  The  VORPAL  simulation  is  a  fully  periodic  box,  containing  atomic 
Oxygen,  0+  and  e'.  The  evolution  of  the  different  charge  states  is  shown  in  Figure  2 
(right).  For  the  given  electron  energy,  the  first  charge  state  not  present  in  the  initial 
distribution  is  going  to  be  02+,  but  as  soon  as  it  is  available  in  high  enough  concentration, 
the  next  ionization  state  is  being  produced  and  the  ionization  rate  of  02+drops.  The 
overall  spectrum  obtained  by  integrating  the  ionization  rates  in  time  is  shown  in  the  left 
panel.  For  comparison,  the  measured  current  at  VENUS  [8]  is  shown.  Good  agreement  is 
obtained  for  5  charge  states.  As  mentioned  previously,  the  charge  state  M/Q=4  is  mixed 
in  the  measured  spectrum.  We  therefore  plotted  a  token  value  of  0  at  M/Q=4.  The  good 
agreement  of  the  simulated  with  the  measured  spectrum  indicates  that  the  kinetic 
ionization  model  in  VORPAL  works  as  expected. 

The  energy  dependent  cross-sections  for  the  ionization  model  in  VORPAL  are  based  on 
the  semi-empirical  model  by  Shull  and  Van  Steenberg  [9],  which  originates  from  the 
astrophysics  community.  This  models  predicts  the  energy  dependent  cross-sections  of 
light  elements,  including,  C,  N,  O,  Ne,  Mg,  Si,  S,  Ar,  Ca,  Fe,  and  Ni. 
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Figure  2:  Left:  Comparison  of  an  oxygen  spectrum  computed  with  the  kinetic  ionization  model  in 
VORPAL  [bars]  with  an  observed  spectrum  [red  stars].  The  measurement  at  Oxygen  4+  was  blended  with 
Helium  and  therefore  not  shown  in  the  picture.  Right:  Temporal  evolution  of  the  different  charge  states  of 
Oxygen  for  a  population  starting  with  pure  Oxygen,  and  0+:  02+  (black),  03+  (brown),  04+  (organge),  05+ 
(yellow),  06+  (purple),  07+  (blue). 


2.2  Fusion  Reaction 


The  first  part  of  this  project  was  to  implement  a  cross-section  model  for  fusion  reactions 
into  VORPAL.  Due  to  the  relatively  large  cross-section  at  low  energies,  we  chose  the 
Deuterium-Tritium  (D-T)  reaction  for  all  further  studies.  We  chose  a  parameterized 
cross-section  model  [6]  of  the  form 


oy(£)  = 


A5  + 

(Aa-AJE?+ 4 

'a2 

E 

exp(^41£,“12)  -1 

} 

(1) 


with  Ai  =  45.95  bam  keV,  A2  =  50200  bam  keV,  A3  =  3.98e-3  keV'1,  A4=1.297  and 
A5=409  bam  keV,  and  the  total  cross  section  aT  measured  in  bam  (10~24  cm2)  and  the 
energy  E  in  keV. 

Figure  3  shows  the  cross-sections  as  a  function  of  energy.  The  peak  cross-section  is  on 
the  order  of  10~27  m2,  centered  around  an  energy  of  130  keV. 

One  of  the  main  differences  between  the  kinetic  ionization  model  and  the  fusion  model  is 
the  magnitude  of  the  cross-section,  which  is  about  7  orders  of  magnitude  smaller  for 
fusion  than  for  ionization.  We  therefore  did  not  track  the  generated  4He  (a  particles)  in 
the  simulation,  but  rather  estimated  the  number  of  neutrons  generated  in  the  interaction 
based  on  the  particle  statistics.  This  de-coupled  approximation  is  valid  because  the  a- 
particles  are  too  few  to  have  an  impact  on  the  plasma  dynamics  (i.e.  energy  deposition), 
despite  the  short  mean  free  path  at  such  high  densities. 
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2.3  Fusion  of  Shock-Accelerated  Ions 

It  has  previously  been  reported  that  the  interaction  of  a  femto-second  laser  pulse  with 
over-dense  plasmas  can  lead  to  the  formation  of  a  shock  wave  in  the  target  and  to  the 
acceleration  of  ions  at  this  shock  front  [10].  In  a  target  consisting  of  a  D-T  mix,  this  will 
therefore  lead  to  a  population  of  D-T  ions  streaming  through  a  D-T  background.  If  the 
energy  of  the  shock-  accelerated  ions  is  large  enough  compared  to  the  background  ions, 
one  expects  that  fusion  reactions  and  therefore  neutron  generation  will  occur  in  the 
vicinity  of  the  shock  front 

In  order  to  test  this  hypothesis,  we  modeled  the  interaction  of  a  femto-second  laser  pulse 
with  a  slightly  over-dense  D-T  target.  We  assumed  a  plasma  of  density  2.8T0  m'  being 
irradiated  by  a  60  fs,  0.8  pm  laser  pulse  with  an  intensity  of  1=1.8  1017  W/cm2  and  a  spot 
radius  of  3  pm.  The  plasma  was  over-dense  with  a  ratio  of  n/ncr  =  1.6.  The  plasma  was 
modeled  in  a  2D  configuration  on  a  3000  x  5  cell  grid.  The  grid  resolution  in  longitudinal 
direction  is  1.6  nm  and  10  nm  in  transverse  direction.  Perfectly  matched  layer  (PML) 
boundary  conditions  were  employed  along  both  directions  of  propagation  and  periodicity 
was  assumed  in  the  transverse  direction. 


D-T  Crossections 

1 5 1  - - -  - T_ 


Figure  3:  Energy  dependent  cross-section  for  the  D-T  fusion  reaction  based  on  the  model  in  [6],  The  peak 
cross-section  is  around  100  keV  and  reaches  10'27  m2  (10  barn). 

Figure  4  shows  phase  space  projections  of  the  interaction  at  different  times.  Due  to  the 
lighter  mass,  the  Deuterium  ions  get  accelerated  to  higher  velocities  than  the  Tritium 
ions.  The  Target  Normal  Sheath  Acceleration  (TNSA)  mechanism  accelerates  both  D  and 
T  ions  at  the  backside  of  the  plasma  (at  x=2.5  pm).  However,  in  addition  to  the  TNSA, 
particles  are  accelerated  in  the  interior  of  the  plasma  at  the  shock  front.  This  results  in  a 
fast  population  of  D  (and  T)  propagating  forward  through  the  background  plasma. 
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Figure  4:  Phase  space  projections  of  the  D  (blue)  and  T  (red)  ions  at  different  times:  (a)-  96fs,  (b)-0.48ps 
and  (c)-  0.96ps. 
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Figure  5  shows  the  energy  distribution  of  the  ions  for  the  above  parameters.  While  these 
simulations  demonstrate  that  a  fast  population  of  D  can  be  generated  by  the  laser-plasma 
interaction,  the  resulting  ion  energies  are  only  on  the  order  of  50  keV,  for  which  the 
cross-section  (see  Figure  3)  is  too  small  to  generate  any  significant  number  of  fusion 
events. 


100 

80 


Figure  5:  Spatial  energy  distribution  of  the  D  (blue)  and  T  (red)  population  after  0.48ps  (left)  and  0.96ps 
(right)  for  a  laser  intensity  of  1=1.8  1017  W/cm2. 


Figure  6  shows  the  ion  phase  space  and  energy  of  an  interaction  with  increased  laser 
intensity  of  I  =  2.6  10  W/cm  .  All  other  parameters  are  identical  to  the  previous 
simulations.  The  energies  of  the  shock-accelerated  ions  are  well  above  the  120  keV 
required  for  a  significant  number  of  fusion  events  to  occur. 


Figure  6:  Phase  space  projection  (left)  and  ion  energy  (right)  for  the  interaction  of  a  60  fs,  2.6  1017  W/cm2 
laser  pulse  with  a  D-T  plasma. 


Using  the  newly  developed  fusion  model  in  VORPAL,  we  were  able  to  investigate  the 
neutrons  generated  in  this  interaction.  Figure  7  shows  the  temporal  evolution  of  these 
neutrons.  The  peak  neutron  flux  occurs  when  the  shock-accelerated  ions  propagate 
through  the  unperturbed  plasma.  Once  they  reach  the  former  back  side  of  the  plasma, 
which  is  now  accelerated  via  TNSA,  the  relative  velocity  is  reduced  and  the  number  of 
fusion  events  is  therefore  reduced. 
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Figure  7:  Left:  Temporal  evolution  of  the  neutron  flux  generated  in  the  interaction  of  a  60  fs  laser  pulse 
with  a  D-T  target.  Right:  Scaling  of  the  peak  neutron  flux  with  increasing  laser  intensity. 


The  energy  of  the  shock-accelerated  ions  is  highly  dependent  on  the  shock  velocity, 
which  is  in  turn  dependent  on  the  laser  intensity.  As  these  ions  are  responsible  for  the 
neutron  generation,  the  laser  intensity  will  control  the  neutron  flux.  Figure  7  shows  the 
scaling  of  the  peak  neutron  flux  as  a  function  of  laser  intensity.  At  intensities  >  3  1017 
W/cm2,  the  target  was  destroyed. 

In  summary,  using  the  newly  developed  fusion  model  in  VORPAL,  we  demonstrated  that 
the  interaction  of  high-intensity  laser  pulse  with  an  over-dense  D-T  target  can  lead  to  the 
generation  of  neutrons. 
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Figure  8:  Reverse  acceleration  of  ions  in  a  two  material  configuration.  Laser  pulse  entered  from  the  left 
and  accelerated  the  first  population.  At  the  interface  of  the  two  materials,  plasma  gets  accelerated  in 
backward  direction. 
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2.4  Reverse  Acceleration  of  Ions 

While  we  demonstrated  in  the  previous  paragraph  that  the  neutron  flux  generated  by  the 
interaction  of  a  femto-second  laser  pulse  with  a  D-T  mix  is  sufficient  for  radiographic 
applications,  one  of  the  drawbacks  of  that  configuration  is  the  significant  fraction  of  the 
D-T  mix  accelerated  in  the  reverse  direction  (i.e.  opposite  to  the  initial  laser  propagation 
direction)  which  is  lost  to  neutron  generation.  These  back-propagating  particles  can  be 
seen  in  the  phase-space  plots  of  Figure  4b-c,  for  x  <  1 .5  pm.  A  more  efficient  target 
design  would  take  advantage  of  that  ablated  material  as  well. 

One  way  of  taking  advantage  of  the  ablated  material  is  to  optimize  the  target  such  that 
ions  are  only  accelerated  in  the  reverse  direction.  This  can  be  accomplished  by  a  target 
consisting  of  an  under-dense  coating  applied  to  an  over-dense  material. 

Figure  8  shows  the  ion  phase  space  of  such  a  two-component  target,  demonstrating  the 
ion  acceleration  mainly  in  the  reverse  direction.  This  is  an  interesting  and  potentially 
important  mechanism  for  this  type  of  target,  which  was  not  previously  investigated. 

This  mechanism  is  not  limited  to  pseudo  ID  simulations,  but  also  works  in  2D.  Figure  9 
shows  the  propagation  of  a  1021  W/m2,  60  fs,  0.8  pm  laser  pulse  in  a  Helium  coated 
aluminum  target.  The  simulation  was  run  on  a  640  x  320  cell  grid,  spanning  a  64  x  32  pm 
domain.  The  neutral  Helium  is  modeled  as  a  gas  at  density  of  1026  #/m'3.  The  propagation 
of  the  laser  pulse  ionizes  the  coating  and  eventually  accelerates  the  ions.  This 
demonstrates  that  the  reverse  acceleration  of  ions  also  works  in  higher  dimensions. 


Figure  9:  Ion  phase  space,  x-y-vx,  showing  an  iso-density  surface  for  He2+  ions  (yellow)  and  Al+  ions 
(green).  Negative  ion  velocities  are  plotted  upwards.  The  laser  pulse  initially  propagates  in  positive  x 
direction,  gets  reflected  off  the  A1  mirror  and  leaves  the  domain  in  negative  x  direction.  The  He2+  ions 
initially  coat  the  inside  of  the  spherical  A1  mirror  and  are  accelerated  in  reverse  direction  by  the  laser  pulse. 


8 


Tech-X  Corporation 


Laser-Fusion  Final  Report 


2.5  Focusing  of  Reverse  Accelerated  Ions 

In  this  section  we  therefore  investigate  whether  the  material  ablated  in  the  interaction  of  a 
strong  laser  pulse  with  an  over-dense  target  can  lead  to  a  neutron  flux  sufficiently  large 
for  radiographic  applications.  An  alternative  application  for  these  reverse-accelerated 
ions  would  be  to  use  them  as  a  driver  for  e.g.  Inertial  Confinement  Fusion.  In  particular, 
one  may  envision  a  combination  of  forward-accelerated  (sheath  mechanism)  and  reverse- 
accelerated  ion  beams  created  by  the  same  laser  pulse,  operating  on  different  sides  of  a 
target,  to  obtain  synchronized  beam-beam  interactions.  Additionally,  one  may  also 
envision  the  target  as  spherically  shaped,  resulting  in  converging  beams  to  increase 
density  and  temperature  at  the  beam  focal  point;  the  target  would  then  include  holes 
through  which  the  laser  beam  can  pass,  similar  to  Hohlraum  targets. 

We  therefore  start  with  investigating  the  possibility  of  converging  ion  beams  from  shaped 
targets.  We  ran  simulations  of  a  10  pm  thick  foil,  with  a  "dimple"  on  the  near  side  having 
a  radius  of  curvature  of  5  pm.  The  plate  is  comprised  of  a  neutral  plasma  with  equal 
numbers  of  electrons  and  singly-ionized  deuterium  ions.  Both  plasmas  have  a  density  of 
2.5  10  m'  .  The  plasma  frequency  is  <x>p=2.825T0  rad/s  and  the  initial  temperature  is 
130’000  K,  resulting  in  a  Debye  length  of  0.866  nm.  While  this  leads  to  numerical 
heating  of  the  plasma,  we  anticipated  that  the  violent  laser-plasma  interaction  will 
dominate  and  provide  at  least  some  qualitative  results.  For  more  quantitative  results,  we 
then  performed  subsequent  simulations  at  higher  resolution. 
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Figure  10:  Flow  field  of  reverse  accelerated  ions. 


The  laser  has  a  wavelength  of  0.8  pm,  for  which  the  plasma  is  slightly  over-dense, 
n/ncnt=1.44.  The  laser  field  has  a  peak  amplitude  of  1.5  1012  V/m  (corresponding  to  an 
intensity  of  3 TO17  W/cm2)  and  a  pulse  duration  of  20fs. 

These  initial  simulations  were  run  on  a  3000x2000  cell  grid  spanning  a  physical  domain 
of  30  pm  x  20  pm  (dx  =  dy  =  10  nm),  and  a  time-step  of  1.73  10~17s.  Perfectly  matched 
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layer  boundary  conditions  were  used  along  the  laser  propagation  direction  and  periodic 
boundaries  in  transverse  direction. 

Figure  10  shows  the  Deuterium  flow  field  after  1.73  ps  (lOO’OOO  timesteps).  The  figure 
clearly  shows  the  backward  motion,  as  well  as  a  focusing  flow  of  the  ablated  material, 
which  could  eventually  lead  to  a  high  density  population  in  the  center  of  the  ‘dimple’. 

One  of  the  key  problems  encountered  in  these  simulations  is  the  long  timescale  on  which 
the  interaction  takes  place.  In  addition,  it  would  be  highly  useful  to  be  able  to  fully 
resolve  high  density  plasmas  to  avoid  numerical  heating,  a  common  problem  in  Particle- 
In-Cell  (PIC)  simulations  when  the  Debye  length  is  not  well  resolved.  In  order  to  further 
investigate  the  focusing  and  to  actually  observe  the  formation  of  the  focused  ion 
population,  we  ran  additional  simulations,  now  with  a  finer  grid  resolution  of  6000x4000 
cells  (dx  =dy  =  5nm)  and  a  time-step  of  1.17'10'17s.  In  addition  to  the  finer  grid 
resolution,  we  also  applied  a  [1  2  1]  box-car  filter  to  smooth  each  component  of  the 
current  deposited  by  the  charges  at  every  time  step.  This  resulted  in  a  significant 
reduction  of  numerical  heating  and  therefore  enabled  simulations  on  long  timescales. 
Finally,  we  increased  the  laser  pulse  duration  from  20  fs  to  50  fs. 


Figure  11:  Two  different  views  of  the  ion  density  after  1.18  ps.  The  laser  pulse  entered  from  the  left,  was 
reflected  at  the  spherical  mirror  and  left  the  domain  to  the  left.  The  reversely  accelerated  ions  form  a  high 
density  population  in  the  center  of  the  mirror.  The  density  of  this  population  is  significantly  higher  than  the 
TNSA  accelerated  population  on  the  right  hand  side  of  the  foil.  The  boundaries  of  the  initial  distribution  are 
shown  as  dotted  lines. 


The  basic  steps  of  the  laser-matter  interaction  in  these  fully  2D  simulations  are  essentially 
the  same  as  in  the  pseudo  ID  simulation  in  Section  2.3;  i.e.,  TNSA.  In  this  mechanism, 
the  electric  field  of  the  laser  pulse  accelerates  the  electrons  in  the  interaction  region, 
which  are  rapidly  accelerated  into  the  target  by  the  Lorentz  force.  The  electrons  stream 
through  the  foil  and  form  a  high-temperature,  charged  cloud  out  on  the  back  side.  The 
ambipolar  field  created  by  this  cloud  (sheath)  accelerates  the  ions  forward,  while  colder 
electrons  are  pulled  back  into  the  material  (return  current).  If  the  laser  pulse  is  long 
enough,  the  electrons  return  into  the  interaction  region  while  the  laser  pulse  is  still  on, 
being  injected  again  into  the  foil  (electron  recirculation).  At  the  same  time,  the  ions  in  the 
interaction  region  start  to  move  in  reverse  direction,  due  to  the  local  positive  charge  left 
inside  the  target  (reverse  accelerated  ions).  See  Figures  1 1  and  12. 
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While  ID  simulations  capture  all  these  effects,  they  are  unable  to  answer  questions  about 
the  spatial  distribution  of  the  reversely  accelerated  ions,  especially  in  a  shaped  target. 
They  also  do  not  properly  capture  effects  from  transverse  gradients;  for  example,  the 
return  current  at  the  periphery  of  the  channel  (outside  the  laser  pulse  envelope)  can 
generate  strong  magnetic  fields  and  affect  the  ion  propagation. 


Figure  12:  Evolution  of  the  ion  density  over  time.  The  laser  pulse  propagated  along  the  x  direction  and  got 
reflected  at  the  spherical  mirror.  Time  evolves  in  negative  Y  direction.  The  ‘ripples’  on  the  ‘canyon-wall’ 
show  the  strong  filamentation  of  the  ablated  material.  After  about  1  ps,  a  high  density  population  in  the 
center  of  the  mirror  has  formed,  propagating  in  backward  direction.  This  population  shows  strong 
divergence,  most  likely  due  to  the  high  local  space  charge. 


2.6  D-T  Acceleration  in  Shaped  Targets 

While  strong  filamentation  of  the  reverse-accelerated  ion  beam  may  prevent  a  good 
focusing  for  beam-beam  interactions,  we  still  need  to  examine  the  potential  use  of  the 
reverse-accelerated  material  shown  in  the  previous  section  as  a  driver  for  nuclear  fusion. 
We  tested  a  generic  configuration  where  the  reverse-accelerated  ions  can  be  focused  onto 
a  target.  Figure  13  shows  the  overall  setup:  the  main  target  consists  of  a  10  pm  thick  tin 
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foil  at  a  density  of  2.510  m'  with  a  5  pm  hemi-spherical  cavity  in  the  center,  which 
acts  as  a  mirror.  The  cavity  is  coated  with  an  under-dense,  2.5  pm  thick  Deuterium  layer 
at  a  density  of  4.4  1026  m"3.  Suspended  at  the  focus  of  the  mirror  is  a  Tritium  ball  with  2 
pm  diameter,  also  at  a  density  of  4.410  m'  . 

The  0.8  pm  laser  has  a  spot  radius  of  2.5  pm  and  an  intensity  ranging  between  1.5  1021 
Wm'  and  2.4  10  Wra'  corresponding  to  peak  electric  fields  of  1.5  TV/m  to  6  TV/m. 
The  simulation  domain  is  30  x  20  pm  covered  by  a  computational  grid  of  6000  x  4000 
cells  and  a  simulation  time-step  of  1.0610~17s.  Perfectly  matched  layer  boundary 
conditions  are  applied  in  the  laser  propagation  direction  and  periodic  boundary  conditions 
in  the  transverse  direction.  Similar  to  the  simulations  in  the  previous  section,  we  chose  an 
electron  temperature  which  results  in  a  Debye  length  resolved  by  the  computational  grid, 
as  well  as  filtering  of  the  currents  in  order  to  suppress  numerical  noise. 


Figure  13:  Setup  of  shaped  target  configuration:  A  high  density  shaped  tin  target  (red)  is  coated  with  an 
under-dense  deuterium  film  (green)  with  a  Tritium  ball  (blue)  suspended  in  the  center  of  spherical  mirror. 
The  Tritium  ball  has  a  diameter  of  2  pm,  the  Deuterium  coating  has  a  thickness  of  2,5  pm.  The 
cavity/mirror  radius  is  5  pm.  The  laser  pulse  enters  from  the  left,  is  reflected  off  the  spherical  mirror  and 
gets  focused  in  the  center  of  the  mirror. 


Figure  14  shows  the  development  of  the  high  density  jet  of  ablated  Deuterium.  The  over- 
dense  lead  foil,  coated  with  the  deuterium,  is  not  shown.  The  peak  density  in  the  ejected 
Deuterium  can  reach  up  to  1.5  times  the  original  Deuterium  density.  The  Tritium  ball 
explodes  as  a  result  of  expulsion  of  electrons  by  the  laser  pulse  (Coulomb  explosion)  and 
interacts  with  the  counter-streaming  Deuterium. 
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Figure  14:  Time  evolution  of  the  D  and  T  density  for  a  laser  amplitude  of  E=3el2V/m  (left)  and  6el2  V/m 
(right)  after  O.lps  (a),  0.3ps  (b),  0.6ps  (c)  and  0.9ps  (d). 
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Figure  15:  Deuterium  (left)  and  Tritium  (right)  density  distribution  after  1.33  ps  for  a  peak  electric  field  of 
E=  3  TV/m. 


Figure  15  shows  the  Deuterium  and  Tritium  density  distribution  after  1.33  ps.  While  at 
that  time  the  reverse  accelerated  ions  form  a  high  density  peak  at  the  original  location  of 
the  Tritium  ball,  the  Tritium  has  been  spread  over  a  large  volume  due  to  Coulomb 
explosion. 


2.7  Neutron  Generation  in  Shaped  Target 

While  the  previous  paragraph  has  shown  that  the  ablated  material  in  a  shaped  target  can 
be  focused  and  generate  a  high  density  peak,  it  remains  to  be  seen  whether  the 
accelerated  ions  have  the  right  energy  to  undergo  fusion. 

Figure  16  shows  the  temporal  development  of  the  neutron  generation  inferred  from  the 
Deuterium  and  Tritium  densities  and  local  mean  velocities  for  two  different  laser 
amplitudes.  For  both  amplitudes,  neutron  generation  initially  happens  where  the  ions 
from  the  Coulomb  explosion  of  the  Tritium  ball  penetrate  the  Deuterium.  Later,  the 
interaction  of  the  Deuterium  jets  creates  another  location  of  high  neutron  yield.  Finally, 
the  Tritium  density  at  the  original  Tritium  ball  location  becomes  so  low,  that  neutrons  are 
no  longer  produced  there. 

Integrating  the  number  of  generated  neutrons  over  the  entire  simulation  domain  allows 
estimating  the  total  neutron  flux  as  a  function  of  time.  Figure  17  shows  the  temporal 
evolution  of  the  neutron  flux  for  different  laser  peak  electric  fields.  While  increasing  the 
laser  peak  field  leads  quicker  to  the  generation  of  neutrons,  the  ions  exceed  the  drift 
velocity  of  peak  fusion  cross-section.  They  therefore  are  lost  to  the  generation  of 
neutrons.  On  the  other  hand,  for  too  weak  a  laser  field,  the  accelerated  ions  are  not  fast 
enough  and  therefore  do  not  lead  to  nuclear  fusion.  For  the  given  configuration,  a  peak 
electric  field  of  3  TV/m  maximizes  the  neutron  flux. 
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Figure  16:  Time  evolution  of  the  neutron  generation  for  a  laser  amplitude  of  E=3el2V/m  (left)  and  6el2 
V/m  (right)  after  0.3ps  (a),  0.6ps  (b)  and  0.9ps  (c). 
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Figure  17:  Neutron  flux  evolution  for  laser  amplitudes  E=1.5  TV/m  (dash-dotted),  E=2  TV/m  (dash-triple- 
dotted),  E  =  3  TV/m  (solid)  and  E  =  6  TV/m  (dashed).  Increasing  the  amplitude  leads  to  an  increased 
momentum  of  the  accelerated  ions  and  therefore  to  a  relative  ion  energy  which  is  beyond  the  peak  cross- 
section  for  the  D-T  fusion  reaction. 


The  neutron  flux  in  Figure  17  is  computed  at  a  distance  of  1  m  away  from  the  target. 
Even  assuming  a  reduction  of  the  flux  by  about  two  orders  of  magnitude  due  to  an 
external  moderator  should  lead  to  a  sufficient  flux  of  neutrons  for  radiographic 
applications. 


3  Summary 

In  summary,  we  found  that  ions  ablated  in  the  interaction  of  a  strong  laser  pulse  with  an 
over-dense  plasma  can  be  focused  by  shaping  the  target.  In  case  of  a  spherically  shaped 
target,  this  leads  to  the  formation  of  a  high  density  ion  population  in  the  center  of  the 
spherical  mirror.  Target  shaping  has  been  successfully  used  to  focus  ions  accelerated  via 
TNSA  on  the  backside  of  the  target.  These  accelerated  ions  are  then  used  as  drivers  for 
Inertial  Confinement  Fusion.  Our  simulations  indicate  that  a  similar  population  can  be 
generated  at  the  interaction  site,  leading  to  a  drive  beam  traveling  in  the  backward 
direction. 

We  have  shown  that  the  interaction  of  a  femto-second  laser  pulse  with  a  slightly  over- 
dense  target  consisting  of  a  D-T  mix  can  be  used  to  generate  a  neutron  flux  sufficiently 
large  for  radiographic  applications.  However,  a  large  fraction  of  the  D-T  mix  does  not 
participate  in  the  D-T  fusion  reaction  and  is  therefore  lost. 

In  a  last  configuration,  the  laser  pulse  was  used  to  strip  a  Tritium  ball  from  all  its 
electrons.  The  resulting  Coulomb  explosion  was  used  to  accelerate  the  Tritium  ions.  The 


16 


Tech-X  Corporation 


Laser-Fusion  Final  Report 


reverse  ion  acceleration  mechanism  in  a  shaped  target  generates  a  counter-streaming 
population  of  Deuterium  ions.  Depending  on  the  peak  electric  field  of  the  laser  pulse,  a 
neutron  flux  large  enough  for  radiographic  applications  can  be  generated. 

The  above  simulations  pushed  40  million  particles  on  a  6000x4000  cell  grid  for  about 
110,000  time-steps  and  produced  about  250  GB  of  data  and  took  about  500  CPU  days 
(12,000  CPU  hours).  Simulations  of  this  size  are  too  large  for  performing  e.g.  Parameter 
studies.  One  possibility  to  reduce  the  simulation  time  is  to  employ  higher-order  particle 
shapes,  which  interpolate  the  electromagnetic  field  not  only  from  the  surrounding  cell, 
but  from  a  more  extended  region  in  space.  Using  the  field  values  from  three  neighboring 
cells,  it  is  possible  to  relax  the  constraint  to  resolve  the  Debye  length  and  only  requires 
resolution  of  the  electron  inertial  length.  This  results  in  a  grid  resolution  which  is  lower 
by  a  factor  (c/vtherm)3  which  can  easily  exceed  a  factor  of  1,000.  In  addition,  this  yields  to 
an  increased  maximum  possible  time  step  due  to  its  coupling  via  the  CFL  condition.  Time 
savings  of  factors  10,000  and  more  can  therefore  be  achieved. 


4  Conclusions 

In  this  project,  we  have  implemented  a  kinetic  model  for  nuclear  fusion  in  the  plasma 
simulation  code  VORPAL.  We  therefore  extended  a  kinetic  reaction  model  implemented 
in  a  different  project  by  cross-sections  for  the  Deuterium-Tritium  nuclear  reaction. 

We  investigated  the  interaction  of  a  medium  powered  laser  pulse  with  a  slightly  over- 
dense  thin  Deuterium/Tritium  target.  The  interaction  leads  to  the  formation  of  a  shock 
wave  inside  the  plasma,  which  can  accelerate  particles  at  the  shock  front.  This  leads  to 
the  formation  of  a  population  of  D  (and  T)  streaming  relative  to  a  D/T  background. 

Using  our  newly  developed  fusion  reaction  model,  we  were  able  to  compute  the  neutron 
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flux  generated  in  this  interaction.  The  peak  neutron  flux  was  on  the  order  of  10  #/cm  s, 
with  a  strong  dependence  on  the  incident  laser  pulse  intensity. 

While  the  overall  number  of  particles  per  shot  is  relatively  small  due  to  the  short 
interaction  time,  the  high  repetition  rate  of  commercially  available  lasers  with  the 
required  parameters  can  lead  to  high  average  neutron  numbers. 

Discussions  with  experimentalists  have  indicated  that  slightly  over-dense  thin  targets  can 
be  realized  using  foams  or  gels  or  even  using  high-density  gas  jets.  The  simulation  results 
therefore  could  be  tested  experimentally  and  possibly  lead  to  the  development  of  a 
compact  neutron  source. 

We  also  performed  simulations  of  reverse  acceleration  by  laser-plasma  interaction  in 
shaped  targets.  We  have  demonstrated  that  a  focusing  of  the  ablated  material  occurs  and 
leads  to  a  high-density  ion  population  propagating  in  reverse  direction.  This  population 
could  be  used  as  a  driver  for  Inertial  Confinement  Fusion. 


17 


Tech-X  Corporation 


Laser-Fusion  Final  Report 


Particle-in-Cell  simulations  of  high-density  plasmas  are  computationally  highly 
demanding  and  algorithmic  enhancements  should  be  investigated.  For  example,  higher 
order  particle  shapes  have  demonstrated  that  considerable  time  savings  can  be  achieved 
compared  to  first  order  interpolation  for  particles. 
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